SO(5) Symmetry in t—J and Hubbard Models 
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Numerical and analytical results are reviewed, which support SO(5) symmetry as a concept uni- 
fying superconductivity and antiferromagnetism in the high-temperature superconductors. Exact 
cluster diagonalizations verify that the low-energy states of the two-dimensional t—J and Hubbard 
models, widely used microscopic models for the high-T c cuprates, form SO(5) symmetry multiplets. 
Apart from a small standard deviation (~ J/10), these multiplets become degenerate at a critical 
chemical potential /x c (transition into doped system). As a consequence, the d— wave superconduct- 
ing states away from half- filling are obtained from the higher spin states at half-filling through SO (5) 
rotations. Between one and two dimensions, using weak-coupling renormalization, a rather general 
ladder Hamiltonian including next-nearest-neighbor hopping can be shown to flow to an SO(5) sym- 
metric point. Experimental tests and consequences such as the existence of a 7r-Goldstone mode both 
in the insulator and superconductor and, in particular, the relationship between the photoemission 
spectra of the insulator and superconductor, are emphasized. 



I. INTRODUCTION 

As their most prominent universal feature, high- 
temperature superconductors (HTSC) always display an- 
tiferromagnetism (AF) and d— wave superconductivity in 
close proximity, in their phase diagram. Recently, a uni- 
fying theory has been proposed, according to which these 
two at first sight radically different phases are "two faces 
of one and the same coin" . They are unified by a com- 
mon symmetry principle, the SO(5) symmetry [fjj. In 
the meantime, numerical simulations, i.e. exact cluster 
diagonalizations, have verified that this symmetry prin- 
ciple is obeyed for widely used microscopic models for 
the high-T c cuprates, namely the (two-dimensional) t—J 
and Hubbard models, to within an accuracy which is bet- 
ter than the average superconducting energy gap |^| |J. 
These microscopic models are known to reproduce salient 
experimental results of the normal state of the high-T c 
cuprates, including in particular the antiferromagnetic 
state. SO(5) symmetry then implies that these micro- 
scopic models also display d— wave superconductivity. 
This gives rise to a microscopic description of the complex 
phase diagrams of the HTSC from the insulating antifer- 
romagnetic phase over to the metallic normal state and, 
finally, to the d— wave superconducting phase. These re- 
sults were further corroborated between one and two di- 
mensions by weak-coupling renormalization-group calcu- 
lations, which demonstrated that rather general ladder 
Hamiltonians, including longer-ranged interactions m 
and hoppings flow to an SO(5)-invariant fixed point 
0. In this paper, the basic numerical and analytical cal- 
culations are summarized, and recent insights regarding 
the microscopic principle behind SO(5) symmetry J8j are 
reviewed. 

The paper is organized as follows: 

Section [n| describes a new idea and dynamical principle 
behind SO (5) symmetry. It is shown that a formulation 
in terms of triplet and hole fluctuations around an "RVB 



vacuum" allows for a physically transparent demonstra- 
tion of the corner stone in SO(5) theory, i.e. that AF and 
SC are "two faces of one and the same coin" . By start- 
ing from this "RVB vacuum" |f2 >, which represents the 
spin liquid state at half-filling, we demonstrate that an 
AF ordered state can be generated by forming a coherent 
state 



which corresponds to z-like triplets condensed into the 
q—TT, i.e. AF wave- vector state. However, in the SO (5) 
theory, the z-like triplet with momentum n and the hole 
pair with momentum are components of one and the 
same SO(5) vector. They are rotated into each other 
by the SO (5) generating operator tt. This implies that 
the above coherent state with condensed triplets can be 
SO(5)-rotated into a corresponding coherent state with 
t\(q = tt) replaced by the (hole-) pair creation opera- 
tor A^(g = 0). This state corresponds to hole pairs con- 
densed into the q = state, i.e. a superconducting state. 
In other words: both the AF and the SC state can be 
viewed as a kind of condensation out of the RVB state, 
or the spin liquid. If the so-constructed AF state is the 
actual ground state at half-filling, then this physically 
very appealing SO (5) construction yields automatically 
the ground state in the doped situation, i.e. the SC state 



In Section [II 



we summarize numerical evidence for 
the approximate SO (5) symmetry of the two-dimensional 
(2D) Hubbard and t—J model. The SO(5) symmetry or- 
ganizes the low-energy degrees of freedom and gives a 
new microscopic picture of the transition from an AF 
ground state to the d— wave SC state as the chemical po- 
tential is varied Q: Our results show that the d— wave 
SC ground states away from half-filling are obtained from 
the higher spin states at half-filling through SO (5) rota- 
tions. We use a general and direct recipe for checking 
microscopic Hamiltonians for SO (5) symmetry, i.e. the 
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concept of "superspin multiplets" M. The basic idea 
is that the low-energy excited eigenstates of a cluster 
display a definite structure characteristic of a particular 
symmetry, a scheme which has already provided convinc- 
ing evidence for the long-range order of the Heisenberg 
AF on a triangular lattice We have numerically di- 
agonalized the low- lying states of the t—J model near 
half-filling and found that they fit into irreducible repre- 
sentations (irreps) of the SO(5) symmetry group. At a 
critical value fi c of the chemical potential, the superspin 
multiplets are nearly degenerate and, therefore, higher 
spin AF states at half-filling can be freely rotated into 
dSC states away from half-filling. Our overall exact- 
diagonalization results, when further combined with a 
detailed spectroscopy of "SO(5)-allowed" and "SO (5)- 
forbidden" transitions between the superspin multiplets, 
suggest that the low-energy dynamics of the t—J model 
fl can be described by a quantum nonlinear a— model 
with anisotropic couplings, and the transition is that of 
a superspin flop transition [jjj. It is truly remarkable that, 
while the physical properties of AF and dSC states ap- 
pear to be diagonal opposite, and they are characterized 
by very different form of order, there exists, nevertheless, 
a fundamental SO (5) symmetry that unifies them. 

In these two-dimensional (2D) microscopic models, the 
undoped situation - in agreement with the experimental 
situation in the cuprates - corresponds to that of a Mott 
insulator with broken SO (3) or spin rational symmetry: 
long-range AF order is realized. The SO(5)-symmetry 
principle then tells us how this long-range magnetic order 
and the accompanying low-energy spin excitations are 
mapped onto the corresponding off-diagonal long-range 
SC order and the low-energy "Goldstone bosons" (the 
7T— mode) in the doped situation 0-^]. 

However, there exists also a second class of Mott-type 
insulators without long-range AF order, i.e. spin liquids, 
which have a gap to spin excitation. There is growing 
experimental evidence that they are also intimately re- 
lated to the physics of high-T c compounds: Not only do 
these compounds show above the Neel temperature and 
superconducting transition temperature at small dopings 
signs of such a spin gap, but there exist also copper-oxides 
with a Cu02 plane containing line-defects, which result in 
ladder- like arrangements of Cu— atoms (for a summary, 
see M ) ■ These systems can be described in terms of cou- 
pled two-leg ladders (l(J , which exhibit a spin gap in the 
insulating compound and thus belong to the spin-liquid 
Mott-insulator variety. Also the related "stripe phases" 
of the 2D Cu02 planes in the cuprate superconductors 
have recently received considerable attention 11 12]. In 



mode of the Mott insulator evolves continuously into the 
"7r— resonance" mode of the superconductor. This SO (5) 
symmetric model offers a reference point around which 
departures from the SO (5) symmetry can be studied. 
Then, two key questions arise, the first being the rela- 
tionship of the exact SO(5) ladder to the "physical" t—J 
or Hubbard ladders and the second regarding the con- 
nection to the other variety of Mott insulators, i.e. the 
ones with long-range AF order. 

With regard to the first question, progress was recently 
made in the regime of weak-coupling: Using the weak- 
coupling renormalization group method, two indepen- 
dent works J||| have demonstrated that rather generic 
ladder models at half- filling flow to an SO (5) symmetric 
fixed point. This work is reviewed in section |y|. 

In section [V, we also summarize work, which tries 



these systems, the apparent connection between the spin 
gap and superconductivity must be explained. 

In order to illustrate how the SO (5) theory can, in 
principle, cope with the challenge, an exactly SO (5) in- 
variant ladder model has recently been constructed [Q. 
In particular, it was shown that the spin-gap magnon 



to numerically attack both questions in the experimen- 
tally relevant intermediate to strong-coupling regimes. 
As shown in this recent work j8j, SO (5) symmetry has 
profound implications for the dynamical correlation func- 
tions, most notably the single-particle spectrum. Specific 
predictions of SO (5), like a "generalized rigid band be- 
havior" (l^Jl^l and the appearance of sidebands in the in- 
verse photoemission spectrum |l5| may indeed have been 
observed long ago in the actual 2D t—J model and also in 
recent angular-resolved photoemission experiments [fl6|| . 
Motivated by the present ladder theory, we have carried 
out more detailed spectroscopies on the 2D model and 
obtained results in strong support of SO (5). 

Finally, we and other authors |j|||||29 have demon- 
strated that, despite the, at first sight, rather unphysical 
parameter values of the SO (5) symmetric ladder model, 
a "Landau mapping" to the more realistic t—J model 
is feasible. This may suggest that the SO(5) symmet- 
ric ladder is indeed the generic effective Hamiltonian for 
2-leg ladder systems and for the above spin-liquid Mott- 
insulator variety, in general. 



II. ONE-TO-ONE CORRESPONDENCE OF 
ANTIFERROMAGNETISM AND 
SUPERCONDUCTIVITY 

In the HTSC, the dominant charge-carrier dynamics 
takes place in the two-dimensional (2D) Cu02-planes 
p6fl . Each Cu02 unit cell contains an effective magnetic 
moment of spin | , essentially due to the Cu ion. Neigh- 
boring Cu-spins form singlets - the energy win due to the 
singlet formation, the magnetic exchange J is relatively 
large ~ 120meV ~ 1400K. On the other hand, the tem- 
peratures for the transition into both low-temperature 
phases, the AF and the SC phases, T^eei and T c , are both 
significantly lower and of similar magnitude (~ 250K for 
T/veei and ~ 100K for T c ). Already this order of magni- 
tude suggests that the mechanism of superconductivity 
does not directly result from the singlet formation, but 
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is instead related to the mechanism, which results in AF 
in the undoped, insulating situation. 

Indeed, both low temperature phases of the cuprates 
are "ordered" : In the undoped case, i.e. in the insulator, 
we have AF order, in the doped case, the phase coherence 
of the superconductor and both phases are in the HTSC 
in immediate vicinity (the "spin glass phase" , which oc- 
curs in some cuprates, is likely due to disorder). There- 
fore, it seems tempting to unify these low-temperature 
phases, despite the fact that on first glance, they appear 
dramatically different: On the one hand, the insulator 
and, on the other hand, the ideal conductor, i.e. the su- 
perconductor. 

Let us consider, at first, the insulator: At high- 
temperatures ~ 1000 K, the singlet pairs are completely 
disordered. This state is termed, therefore, a spin liquid. 
How does one arrive from this disordered state at high 
temperatures, at an ordered Neel state at low tempera- 
tures? To solve this problem, we shall consider a dynami- 
cal principle, which gives a particularly simple and trans- 
parent demonstration of the key feature of SO(5) theory, 
namely the one-to-one correspondence of antiferromag- 
netism and superconductivity: According to this dynam- 
ical principle 0] , the ordered AF state can be considered 
as a kind of Bose-Einstein condensation of triplet exci- 
tations. The SC state, on the other hand, corresponds 
to a Bose-Einstein type of condensation of Cooper pairs, 
namely in the cuprate materials of hole pairs. Triplet ex- 
citation and hole-pairs excitation are "two phases of one 
and the same coin" in similarity to other unifying con- 
cepts such as the isospin theory of proton and neutron 
in nuclear physics p"8[ |. The condensation energy yields 
then a new lower temperature scale T^eei and T c . 

Here, SO (5) symmetry enters in detail: Like a "mag- 
nifying glass" , this symmetry principle allows to differ- 
entiate the low-energy physics of the order of T c and, 
in particular, to unify the, in principle, competing AF 
and SC phases (specifically, this symmetry principle "ro- 
tates" triplet excitations into hole-pair excitations and 
vice versa). Unifying principles via symmetry are, of 
course, known from many fields of physics, such as quarks 
in high-energy physics, the predictions of which were in- 
spired by the SU(3) classification of hadronic spectra. 
How does such a unification take place in the HTSC, 
i.e. in solid-state physics? 

The order parameter of the AF is the sublattice mag- 
netisation, a real 3-dimensional vector; if this vector is 
different from zero, we have AF order. Consider now two 
fixed neighboring sites in the 2-dimensional AF, let's say 
in the configuration (Tl)- This fixed spin configuration 
may be viewed as being due to an superposition of the 
singlet (U - IT) with the (S z = 0)-triplet (|| + IT)- 
To create the macroscopic 2-dimensional AF in the Cu02 
plane, therefore, we have to mix triplet excitations al- 
ready at high temperatures into all possible singlet con- 
figurations of the spin liquid. The AF states then re- 



sult as a kind of "condensation" of the triplet excitations 
into the lowest possible energy state || . The three com- 
ponents of the triplet correspond to the three possible 
orientations of the AF and the density of the "condensed 
triplets" corresponds to the magnitude of the sublatice 
magnetisation. This dynamical principle illustrates the 
relation between AF and SC states rather clearly: If, in 
the AF state (we shall see below that it corresponds to 
a coherent state) the triplet excitation operator is re- 
placed by a hole pair creation operator, we obtain a 
coherent state which creates a macroscopic number of 
Cooper pairs, i.e. a SC state. The "rotation" AF^SC, 
therefore, is described by an operator, the it— operator of 
SO (5) theory |], which replaces triplets by hole pairs. 

In the following, we give a particularly simple illustra- 
tion for this key- feature of SO (5) theory, which can be 
worked out for a ladder || , and discuss the equivalence of 
antifcrromagnetism and superconductivity for this exam- 
ple. The ground state of the ladder models is actually a 
resonating valence bond (RVB) type of vacuum without 
AF long-range-order However, for illustrative pur- 
poses, let us now construct an AF ordered state (which 
is in general not an eigenstate of the Hamiltonian) by 
condensing magnons into the RVB ground state. One 
can express the operator of staggered magnetization in 
z-direction as 

M s = 5> i7ro ( P„(T!) - P„(IT) ), (i) 

n 

where e.g. -P n (TI) projects onto states where the n rung 
has the configuration T| (see Fig. |j). It is easy to see that 

(^n(TI)-Pn(IT))4=C 
( Pn(TI) - Pn(IT) = 5 a , z 4, (2) 

where 4(4«) creates a singlet (a-tiplet) on rung n. 
Therefore 

M s = ^[tt(q = 7r)+t z (q = n)]. (3) 
If we now form the coherent state 

|tt A > = J- e ^l(v=*)\n}, (4) 
V 71 

which corresponds to z-like triplets condensed into the 
k—n state, and treat the t as ordinary Bosons, we obtain 

(*a|A/ s |*a) = V2XN. (5) 

This calculation shows that by starting from a spin liq- 
uid, i.e. an RVB vacuum, an antiferromagnetically or- 
dered state with Ms in z-direction can be generated by 
condensing z-like triplet-excitations into the fc=7r state. 
At this point, we can invoke the SO (5) symmetry of the 
model H , which tells us that since the z-like triplet with 
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momentum tt and the hole pair with momentum are two 
different components of a 5- vector, they are dynamically 
indistinguishable. This means that the AF state, with 
condensed triplets, can be SO(5)-rotated into a state with 
condensed hole pairs. It follows that, if the antiferromag- 
netic state were the ground state at half-filling (which is 
the case for 2D materials and physical ladder systems), 
we can replace all z-like triplets by hole pairs with mo- 
mentum and by SO (5) symmetry automatically obtain 
the ground state in the doped case. The latter then con- 
sists of hole-pairs along the rungs condensed into the k=0 
states and thus is necessarily superconducting. In other 
words: both the antiferromagnetic and the superconduct- 
ing state may be viewed as some kind of condensate 'on 
top of the rung-RVB state. SO(5) symmetry then sim- 
ply implies that the condensed objects are combined into 
a single vector, whence the unification of antiferromag- 
netism and superconductivity follows in a most natural 
way. 

The above derivation makes sense only in a strong cou- 
pling limit, where a ground-state description starting out 
from rung-singlets is appropriate. One might expect, 
however, that similar considerations will apply also for 
cases with a weak coupling within the rungs (||| ( see 
also section IV). 

In two dimensions (2D) the above interpretation of 
SO (5) symmetry hinges crucially on two points: first, 
the excitation spectrum of a 2D spin-liquid or RVB 
state must consist of Bosonic triplets (rather than 
e.g. Fermionic 'spinons'), whose condensation into the 
minimum of their dispersion at (tt, tt) leads to antifer- 
romagnetic ordering. This would imply that the anti- 
ferromagnetic state in 2D also could be interpreted as a 
condensate of triplet Bosons. Work along this line is in 
progress and promising . 

Second, these triplet-like Bosons with momentum 
(tt, tt) must be dynamically equivalent to a d x 2__ y 2 Cooper 
pair. Having established this equivalence one could im- 
mediately conclude that the antiferromagnetic phase (viz 
'condensate of triplets') is identical to the <i-wave su- 
perconducting phase (viz 'condensate of hole pairs') in 
the same sense as two nuclei belonging to the same Iso- 
multiplet are 'identical'. Since the 7r-operator in 2D |l]] 
precisely converts a triplet with momentum (-7T, tt) into 
a d- wave hole pair (see section III), the latter require- 
ment is equivalent to [H, tt] = loqtt. This commutation 
relation, which will be discussed in detail in the next 
section, tells us whether the dynamics of charge carriers 
respects SO (5) symmetry, i.e. the Hamiltonian is in ac- 
cord with the SO (5) rotation tt. Ideally, we should then 
have zero for the commutator. This will happen for a 
critical chemical potential /i c , right where the transition 
into the doped system takes place. Further changing /i 
gives rise to a precession frequency, or energy to perform 
the AF — > d— SC rotation ujq ^ 0. The energy shift u>q 
here would correspond to the mass difference of proton 



and neutron in the Isospin algebra. The validity of this 
eigenoperator relation has already been verified numeri- 
cally §. 



III. SO(5) SYMMETRY AND MICROSCOPIC 
MODELS FOR HTSC 

The dynamical principle of the previous section sug- 
gests to unify the triplet excitations (corresponding 
to a vector) and hole-pair excitations A (described by 
(Rc(A)), Im(A)) into a 5— dimensional vector, the so- 
called superspin vector. More precisely, Zhang suggested 
to group the AF order parameter S(Q) with Q = (tt, tt) 
and d-w&ve SC order parameter A into a single five- 
component vector, the superspin 



» B = (At + A ) S(Q),-i(At-A)) 



(6) 



where A = (i/2) ^2 p (cosp x — cos Py)c p a y c^p denotes 
the d x 2_ y 2 superconducting order parameter. S (Q) = 
S P c Q+ P (TC p stands for the AF Neel vector and a a are 
the Pauli spin matrices. The transition from AF to dSC 
is then viewed as a kind of "superspin flop" transition 
as a function of the chemical potential or doping, where 
the direction of the superspin changes abruptly. This 
transition or SO (5) rotation is formally described by the 
so-called tt— operator. 

The "superspin flop" transition is analogous to the 
problem of a spinning top in a uniform gravitational field, 
or a magnetic moment in a uniform magnetic field, as in 
nuclear magnetic resonance (NMR). As we shall see, in 
the presence of the symmetry-breaking field (which, in 
our case is the chemical potential or doping), the order 
parameter is forced to precess, and this explicit break- 
ing of the SO(5) invariance induces rotations between 
AF and dSC states (with frequency uiq) and governs the 
competition between these two phenomena. 

Mathematically, the f}— operator can straightforwardly 
be constructed as a rotation from AF to dSC, it has 
to patch up the differences in the quantum numbers of 
the corresponding order parameters. Since A has spin 
S=Q, whereas S (Q) has S—l (which is obvious from our 
triplet excitation picture in the last section), tt has to 
carry S=l, i.e. it must be a triplet operator. S (Q) has 
no charge, A has charge ±2, therefore, tt must create 
charge ±2. Finally, S (Q) has momentum Q = (tt,tt) , 
A instead has Q = 0, thus, f} must have momentum 
Q. Combination of these requirements fixes the operator 
uniquely up to a form factor, which is given by rf-wave 
symmetry |B8| . This results in the tt operators (we just 
give one of them) 



7T^ = ^""^ (cos p x — cos p y 



,c p+qT c -pT- 



(7) 



Charge conjugation and spin lowering operation gives five 
other similar operators. In real-space representation, the 
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7r (tt) operator does precisely what we discussed in the 
previous section, namely it replaces triplets oriented in x 
and y planar directions by d x 2_ y 2 hole (electron) pairs. 
These tt operators were first introduced by Demler and 
Zhang p(| to explain the resonant neutron scattering 
peaks in the YBCO superconductors. Together with the 
total spin S a , which is the generator of the SO (3) spin 
rotation forming a subgroup of SO (5), and the charge 
operator Q generating the U{\) charge symmetry (which 
is also a subgroup of SO(5)), the six n (n^) operators 
form the ten generators L a b of the SO (5) algebra (for 
definition of L a b, sec Q] ). As constructed, the tt opera- 
tors indeed rotate the AF order parameter into the dSC 
order parameter, i.e. 



(8) 



and vice versa. 

We have, thus, apparently accomplished the task of 
unifying AF with SC: the corresponding order parame- 
ters are grouped into a five dimensional object, and SC 
is "nothing" but AF viewed in some rotated coordinates 
and vice versa! This construction looks a bit similar to 
the unification of E and B by the Lorentz group. But, so 
far, this is only a mathematical construction, we haven't 
asked if "Mother Nature" approves the SO (5) construct 
or not. 

In the high— T c problem, "Mother Nature" is very com- 
plicated, but we can check the SO (5) symmetry within 
some microscopic models, which are known to reproduce 
salient features of the phase diagram of HTSC. Such mi- 
croscopic Hamiltonians, i.e. Hubbard and t— J models, 
successfully model the Mott-Hubbard insulator to metal 
transition which is driven by the Coulomb cor- 

relation U ~ lOeV, i.e. by "high-energy" physics. They 
also model prominent features of the magnetic interac- 
tions on an energy scale of order J ~ O.leV Jl3|- |l5| , po| pi[ . 
However, their low-energy content of order of the average 
SC gap (~ J/5 — J/10) has so far eluded both analyti- 
cal and numerical investigations JjJ. We shall demon- 
strate that SO (5) symmetry overcomes this major obsta- 
cle; it clarifies the role of competing orders and gives a 
microscopic description of the transition from AF to dSC 
ground states, as the chemical potential is varied. 

One can check the SO (5) symmetry by evaluating 
the commutator between the Hamiltonian with the 
7T— operators. In particular, numerical works, sum- 
marized below, show that the tt operators are approxi- 
mate eigenoperators of the Hubbard Hamiltonian, in the 
sense that 



[H,7t, 



(9) 



where the eigen- frequency uiq, which is the energy to per- 
form the AF— >dSC rotation, is of the order of J, and 
proportional to the number of holes. This relation (H) re- 



minds us of the commutation relation between the trans- 
verse spin components S+ and 5_in a magnetic field and 
the Zeeman Hamiltonian, with loq being proportional to 
the B— field. In our case, the SO(5) symmetry is broken 
explicitly by the chemical potential, i.e. luq is propor- 
tional to the hole count or doping. Thus, the pattern of 
explicit symmetry breaking is simple and familiar, and, 
therefore, easy to handle. 

In the following, we employ a kind of "computer spec- 
troscopy" to test whether the dynamics of the charge car- 
riers, i.e. the Hamiltonian, respects SO (5) symmetry and 
(^) is fulfilled. Ever since the early days of quantum dy- 
namics, group theoretical interpretation of spectroscopy 
revealed deep symmetry and profound unity of Nature. 
Atomic spectra can be fitted into irreducible representa- 
tions (irreps) of SO (3), and the regular patterns which 
emerged from this classification offered fundamental un- 
derstanding of the periodic table. After the discovery of a 
large number of hadrons, the "embarrassment of riches" 
was removed by the classification of hadronic spectra into 
irreps of SU(3) and this hidden regularity inspired the 
predictions of quarks, the fundamental building block of 
the universe. In our quest for understanding the fun- 
damental design of Nature, the importance of symmetry 
can never be over-emphasized. 

In our work, we used a different kind of spectroscopy 
and its classification into a different kind of symme- 
try. The spectroscopy is performed on a computer, 
which numerically diagonalizes microscopic Hamiltoni- 
ans, i.e. Hubbard and t— J models widely believed [^0 21 
to model high-T c superconductors. 

As the simplest 2D lattice model for correlated elec- 
trons, the one-band Hubbard model is defined as 



H 



-t 



£(• 



ct C ia + h.C. ) + 



C/]T 7^7^, (10) 



with nearest-neighbor hopping t and Coulomb correlation 
U. On the other hand, the t— J model has the Hamilto- 



H = P 



-t 



/ , \ C l<7 C j" 



(id) 



p. (11) 



As in the Hubbard model, (ij) denotes a summation over 
nearest-neighbors on the 2D square lattice and P projects 
onto the subspace with no doubly occupied sites. The 
latter constraint reflects the strong correlations in the 
U/t — > oo limit of the Hubbard model. 

First numerical evidence for the approximate SO (5) 
symmetry of the Hubbard model came recently from 
exact diagonalizations of small-sized (10 sites) clusters 
H , studying dynamic correlation functions involving the 
AF/dSC rotation fr operator. We observe that (||) is 
nothing but the ladder operator relation familiar from 
standard Quantum Mechanics problems such as the har- 
monic oscillator or the spin rising (<S+) and lowering (SL) 
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operators. If this equation is fulfilled, then the equal level 
distance luq between the "rungs" of the ladder should ap- 
pear - like in NMR or optical spectroscopy - as a sharp 
peak in the n— ir correlation function, well separated from 
a higher-energy incoherent background. This indeed is 
verified in Fig. |]](a),(b), which displays a typical result 
for the itd correlation function in a ten-site Hubbard clus- 
ter with U/t = 8 for dopings of one hole-pair ((n) = 0.8) 
and two hole-pairs ((n) = 0.6) ||]. We note that ujq is 
a small energy, scaling with the hole count away from 
half-filling ((n) = 1) , i.e. uj = J/2(l - (n)) - 2/m. 




2 4 6 8 10 12 14 



FIG. 1. Dynamic correlation functions of the vTO x \/T0 
Hubbard model with U=8t: (a) 7r^(oj)-spectrum at (n)=0.6, 
(b) 7r^(oj)-spectrum at (n)=0.8, (c) tyI (tj)-spectrum at 
(n)=0.8 

The approximate relation @ is highly nontrivial. One 
could ask if a similar relation would exist for a modified 7r 
operator which rotates AF into s— wave SC order param- 
eters. The answer is negative ||||, as seen in Fig. |l|(c) 
which demonstrates that a tt— rotation with s— wave sym- 
metry just generates an incoherent background and no 
sharp "7T— resonance" . Therefore, there is only an ap- 
proximate symmetry between AF and d— wave SC (and 
not s— wave SC) near half filling. 

In the next step, we use a most general and direct 
recipe for checking microscopic Hamiltonians for SO (5) 
symmetry, i.e. the concept of "superspin multiplets" 
We consider the t— J model, which, because of its 
more limited Hilbert space (no double occupancies), al- 
lows the exact diagonalization of larger systems (18, 20 
sites). Since the t—J model explicitly projects out the 
states in the upper Hubbard band, some of the questions 
p5| , p6]| raised recently about the compatibility between 
the Mott Hubbard gap and SO (5) symmetry can also be 
answered explicitly. In particular, if there is an approxi- 
mate SO(5) symmetry of the microscopic model, the low- 



energy states of this model should fall into irreducible 
representations (irreps) of SO (5). In a given quantum 
mechanical system, the direction of the SO (5) superspin 
vector is quantized in a way similar to an ordinary SO (3) 
spin, and the classically intuitive picture of the preces- 
sion of the SO (5) superspin vector under the influence 
of the chemical potential |i] can be identified with the 
equal level-spacing between the members of SO (5) mul- 
tiplets carrying different charge. Therefore, numerically 
identifying the low-lying states of the microscopic model 
with the SO (5) irreps can lead to detailed understanding 
of the one-to-one correspondence and the level crossing 
between the excited states of the AF and the dSC states, 
and thereby lead us to the microscopic mechanism by 
which the AF changes into the dSC state. While finite- 
size calculations cannot generally be used to prove the 
existence of long-range order in infinite systems, the spec- 
troscopic information about the SO(5) symmetry can be 
used as input for the effective field theory |jl|,^,^8) which 
captures the low-energy and long-distance physics of the 
problem. 

Exact diagonalizations (e.g. jlB)) commonly study 
ground-state correlations, but their spatial decay is often 
inconclusive as a test of order due to small system size. 
Yet it is possible that the (excited) eigenstates show a 
well-defined structure characteristic of a particular sym- 
metry; this provided the convincing evidence for long- 
range order in the spin- 1/2 triangular lattice AF 0. In 
our work, we have pursued exactly such a program in 
exact diagonalizations of the t—J model. 

Consider as a simplest example the precession of 
a spin- 1/2 system in a homogeneous B— field in 
z— direction. B breaks the SO (3) spin-rotation symmetry 
and the spin (expectation value) precesses with Larmor 
frequency around the spatially fixed B— direction. This 
can be read off in a spectroscopic experiment from the 
multiplet structure of the possible spin states (| ± \ >) , 
the degeneracy of which is lifted by the B— field (Zeeman 
effect). In this simplest case, the multiplet structure is 
one-dimensional, extending in ^—"direction". 

In our "computer spectroscopy" , we employ a formally 
analogous recipe for checking the microscopic Hamilto- 
nian for SO (5) symmetry: 

In SO (5), the multiplet structure is two-dimensional, 
i.e. spanning both S^-direction and Q-direction |Q. Q = 
L15 is the total charge, and, in our t—J calculation, it 
stands for the number of doped hole pairs, and, thus, for 
the transition from AF to dSC states. Formally, the mul- 
tiplets are constructed by observing that {S z , Q, C} form 
a set of commuting operators with their quantum num- 
bers labeling states of an SO (5) invariant Hamiltonian 
0. S z is the z— component of the total spin (S z = —L23) 
and C the Casimir operator X) Q <&^a&> which is a nat- 
ural generalization of the total spin operator S 2 . Like 
S 2 in the familiar SO (3) spin-rotation symmetry, C fixes 



G 



the level v of an irreducible representation. Instead of 
S(S + 1) in SO(3), it takes the value v{v + i) for an 
SO(3) level v irreps. Fig. || shows the first four {v = 
to v = 3) irreps of SO (5) with the low- lying states of an 
18-site t—J model with J/t = 0.5, which is a typical J 
value ITBI. 



X ® X X 



x x ® 



A levelv irreps of SO(5) 



X X 
® >< 
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x x cfe x x 
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FIG. 2. The upper diagram illustrates general level v irreps 
of SO(5). Every state can be labeled by Q and S z . The max- 
imal charge is Q = ±v. The states labeled by a x form the 
shape of a diamond, while states inside the nested diamonds 
are labeled by o and A. Overlapping states with same Q and 
S z are distinguished by their S quantum numbers. The lower 
diagrams are for v = 1, 2, 3 irreps of SO(5). The figure shows 
the energies of some low energy states for the 18-site cluster 
with J/t = 0.5. The states are grouped into different multi- 
plets and are labeled by the spin, point group symmetry, and 
total momentum. Ai denotes the totally symmetric, B\ the 
d a .2_ y 2-like representation of the Ca v symmetry group. The () 
symbol denotes as yet unidentified members of the respective 
multiplet 

What is the physical meaning of these multiplets? In 
fact, there is a definite physical meaning behind them, 
which is closely related to our coherent-state description 
of the AF-^dSC transition in section ||. The AF state is 
constructed from the linear superposition of the Q = 0, 
S z = states in each of the level v multiplets. At a 
given level v this state contains v magnons or triplet 
excitations, which correspond to the v— th term in the 
power-series expansion of the coherent-state operator, 



i.e. e 



l + A(tt)-=i + ^-(tt)-=2... 

Correspondingly, the dSC state is constructed from a 
linear superposition of the lowest corner states of each 
level v irreps, and may be viewed as stemming from the 



power-series expansion, e 



AA 



1 + AA'' 



4rA^= 2 . 



From Fig. || we note that the low-lying states indeed 
fit into the irreps of SO (5): all the different quantum 
numbers of the states are naturally accounted for by the 
quantum numbers of the superspin. However, most im- 



portantly, the levels with different charge Q are nearly 
equally spaced! This is explicitly indicated by the sym- 
bol AE in Fig. ||. More precisely, the mean-level spacing 
within each multiplet (up to Q = —2) is —2.9886 with a 
standard deviation of 0.0769. This standard deviation is 
much smaller (~ J/8) than the natural energy scale J of 
the i— J model and comparable to or even smaller than 
the average SC gap. If one now adds the chemical po- 
tential term, = — 2/iQ, then at a chemical potential 
fi c comparable to the mean-level spacing, the superspin 
multiplets are nearly degenerate. In other words, (each 
of the coherent-state contributions to) the AF state can 
freely be rotated in (the corresponding contribution to) 
the SC state. 




FIG. 3. Spectral functions with final states in the Q=— 1 
subspace: dynamical spin correlation function for momentum 
transfer Q, calculated for the 1 Bi(0,0) ground state; spec- 
trum of the 7T -operator, calculated for the 1 ^4i(0,0) ground 
state in the Q=—2 sector; spectra of the 7r-operator, calcu- 
lated for the half-filled 1 y4i(0,0) ground state and the lowest 
half-filled 5 Ai(0,0) state. 

In Fig. ||, this is verified in terms of another "diagnostic 
tool" , which is the spectral function 



A(co) 



1 



1 



7T U! — (H — ti re f) — tV^ 



■0\r/,). (12) 



We see that, if we apply as the operator O the op- 
erator to the "dSC" state in the v = 2 irreps (the 
Q = —2,S Z = state 1 Ai(0, 0)), we end up in pre- 
cisely the same final state (to within computer accu- 
racy) as if we apply the tt operator to the "AF" state 
(Q = 0, S z = 0; 5 A 1 (0,0)). In other words, two suc- 
cessive 7T— rotations rotate us from the AF state to the 
d-SC state. The energy required to perform this rota- 
tion, u)q, is again seen to scale with J (in contrast to a 
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false argument in ref. p5fl). 

Finally, we note that the dynamical spin correlation 
function - where we take as the operator O in ([l2]) the 
magnon operator S(Q) - has a peak at precisely the same 
position as the n— resonance. This confirms the earlier 
conjecture by Demler and Zhang pc|| to interpret a well- 
known spin resonance detected in neutron-scattering ex- 
periments in HTSL, as a 'V— excitation" . More gener- 
ally, this fingerprint of SO (5) symmetry can be described 
as follows: Away from /^ Cj the chemical potential, which 
controls the changes in the total charge, breaks SO (5) 
symmetry. As already mentioned, the effect of this chem- 
ical potential is formally analogous to the effect of the 
B— field on a spin. In the presence of this symmetry- 
breaking field, the order parameter is forced to precess 
and this explicit breaking of the SO (5) invariance induces 
transitions between AF and SC states, and governs the 
competition between these two phenomena. The SO (5) 
theory - in good accord with the experimental data ]3l| ] 
both concerning doping and temperature dependence - 
then identifies this precession frequency wo with a reso- 
nance in neutron scattering. 



IV. SO(5) SYMMETRY IN LADDERS 

A. Introduction 

High-T c materials are antiferromagnetic Mott insula- 
tors displaying long-range AF order at half filling. The 
antiferromagnetic phase is rapidly destroyed upon doping 
and is replaced by the superconducting phase. Below op- 
timal doping (the doping with maximum T c ) and above 
the superconducting temperature there are clear experi- 
mental indications for the opening of a spin gap |MB2[ , 
This phase is termed spin-liquid with properties which 
are quite difficult to reproduce on the theoretical level 
for a two-dimensional system. On the other hand, at and 
close to half filling this spin-gap variety of Mott insulators 
is obtained quite well in half- filled ladder systems: Upon 
doping the ladder, the spin gap persists and the system 
exhibits quasi-long-range d-wave superconducting corre- 
lations, which become dominating in some parameter 
range pX|j . In this sense, ladders systems show properties 
very similar to the phase diagram of cuprate materials, 
the main difference being the fact that correlations are 
"quasi-long-ranged" , i.e. they show power-law behavior, 
since they cannot be truly "long-ranged" because of one- 
dimensionality. Moreover, there also exist copper-oxides 
with Cu02 planes containing line defects, which result 
in ladder- like arrangements of Cu-atoms |3j| . These sys- 
tems can be described in terms of coupled two-leg ladders 
exhibiting a spin gap and thus belong to the spin-liquid 
Mott-insulator variety. Also the related "stripe phases" 
of the 2D Cu02 planes in the cuprate superconductors 



pl| , which have recently received considerable attention, 
can be mapped onto ladder systems. 

It is thus interesting to study the occurrence of SO (5) 
symmetry in ladder systems. As we will discuss, there is 
in fact a natural way to construct an SO (5) symmetric 
model for a two-leg ladder, which has only local interac- 
tions on a rung of the ladder 0. We can thus use the 
ladder system as a theoretical laboratory to check some 
ideas of the SO (5) theory. On the other hand, we want to 
study whether generic models, which are not SO (5) in- 
variant at the bare starting (microscopic) level may show 
SO (5) symmetry in their low-energy regime. An appro- 
priate tool to study this low-energy regime, starting from 
a microscopic Hamiltonian, is the renormalization group 
(RG). This method has proven to be particularly suited 
to study systems between one and two dimensions, for 
example ladder systems |34|-|3(|. 

In the following, we consider generalized Hubbard-type 
ladders with Hamiltonians of the form 

H = -t ( C L C J> + h - c -) - c l,* c i±y,v 

(i,j),a »>" 

+U ^2 n i,t n i,l + v ^2 nin i + ~Y ^2 niTli± y 
+ ^^ Si S i±s . (13) 

i 

(see Figure ||). The sum over i in (|T^) runs over the 
sites of a two-chain ladder, as before, covers pairs 

of nearest-neighbor sites on the same chain and i ± y is 
the nearest-neighbor site to i on the other chain. 




FIG. 4. A generalized Hubbard-type ladder model: the 
Hamiltonian consists of electron hopping terms (t), on-site 
and off-site Coulomb interactions (U, V), and spin-spin ex- 
change terms (J) 

B. Exact SO (5) ladder 

In this paragraph we present "numerical spectroscopy 
experiments" for a recently proposed exactly SO (5) 
symmetric ladder model based on a simplified version 
of the more general Hamiltionian (jf^): The non-local 
Coulomb and spin-spin interactions (the V and J terms 
in fll3|)) are restricted to act within the rungs, but not 
along the legs of the ladder, and hopping is only allowed 
between nearest-neighbor sites (with hopping t along 
the legs and t± within the rungs) (see Fig. [|). Setting 
J± = 4(U+V±), the elementary magnetic and "SC" exci- 
tations of a half-filled single rung, namely the formation 
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of a triplet and the formation of an electron (or hole) pair, 
become degenerate in energy. In this case, one obtains 
an exactly SO (5) symmetric Hamiltonian 0. In a recent 
work ||, this ladder model was studied in detail and re- 
lated to the triplet-excitation picture and coherent-state 
description of section O. Furthermore, this work gives 
the most general construction of the SO(5) irreducible 
representations not only for even numbers of holes (or 
electrons) as in section [II, but also for odd numbers, en- 
countered for example in direct or inverse photoemission 
(PE). Exact diagonalization was used to extract the prac- 
tical application of the selection rules for photoemission 
implied by SO (5) symmetry. 



numbers (p, q) : within a given (p, q) multiplet the 
states are characterized both by S z and the charge Q (in- 
stead of only S z in SO (3)). Fig. || and ^ display some of 
the lowest SO(5) multiplets as well as some allowed PE 
transitions between them (note that for p~q, one gets 
the diamond- like structures shown in Fig. |^). 

These theoretical assertions can be checked numeri- 
cally using the Lanczos technique fl5| . The method 
permits to calculate ground state properties like energy 
and spin expectation values, but also Green's functions 
(cf. (^) ), and, in particular, the photoemission and in- 
verse photoemission spectra (here for a spin-down elec- 
tron): 
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FIG. 5. SO(5) multiplets connected by allowed photoemis- 
sion (PES) and inverse photoemission (IPES) transitions. As 
we inject/remove a .[-electron, PES corresponds to an ar- 
row pointing south-east and IPES to one pointing north-west. 
Here, the initial state is a half-filled state with S z =0 in the 
(1,1) multiplet (open circles). SO(5) selection rules only al- 
low transitions to the (2,1) and the (1,0) multiplets (black 
circles), which are superimposed to the initial (1, l)-multiplet. 
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FIG. 6. PES/IPES from a spin-polarized half-filled state 
(Q=0, S z =+1) and from a doped state (Q=—l, S z =0). Both 
states are members of the (1,1) multiplet like the initial state 
in Fig. H Note that there is no allowed photoemission tran- 
sition to the (1,0) multiplet for either of the two states. 

The classification of SO (3) symmetric states into spin 
multiplets, e.g. a triplet with 5=1, is well known. For 
SO(5) symmetry a similar classification holds, but to 
characterize a multiplet one needs two integer quantum 



^(k,.) = ^(Olct i a;+if _ 1 eo _, o+ c ki 10 ) , 
A IPES (K^) = 1 -^\ Ckl Lo + H _\ o _ tQ+ C t l \0). 

|0) is an initial energy-eigenstate found by Lanczos- 
diagonalization and eo its energy. 

With this "computer experiment" one can study the 
direct and inverse PE within a multiplet structure and 
observe crucial selection rules. The SO (5) multiplets are 
easily identified, because the energies of states belonging 
to one multiplet are degenerate (i.e. identical to com- 
puter accuracy of ca. ICE 14 ). This allows to study the 
evolution of the single-particle spectral function as we 
pass from one multiplet (p, q) to the other, and within 
one multiplet through the different doping and spin lev- 
e\s(Q,S z ). 

The dotted line in Fig. ^ shows the single-particle 
spectrum for the half-filled ground state 1 (0,0)o (where 
the group theoretical notation is a short hand for 
de s enerac y(= s \momentum) Q ); this is an RVB-state of 
rung-singlets (see section |) which actually forms a one- 
dimensional (0, 0) multiplet. Final states can only be- 
long to the 4-dimensional (1,0) irreps ||]. Despite the 
fact that we are using very strong interaction parame- 
ters, there is just one single electron-like band in PES, 
whose cosine-type-of dispersion closely follows the disper- 
sion of a non-interacting electron. The center of gravity 
of this band is given by the energy difference between a 
rung-singlet and a singly occupied rung ||]. 

The solid line in Fig. ^ represents the PES for the half- 
filled 3 (7r, 7r)o state (with S Z ~Q), which carries one triplet 
excitation (magnon) and which belongs to the fivefold 
degenerate (p=l,q=l) multiplet, i.e. the open circles in 
Fig. M. From Fig. || we expect in this case final states 
belonging to both the (2, 1) and the (1, 0) irreps. Indeed, 
in addition to the band seen in the ground state spec- 
tra, which remains practically unchanged, there appears 
a "sideband" close to [i both in PES an in IPES; thus 
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FIG. 7. PES/IPES spectra of a 6- rung ladder (re- 
moval/injection of a [-electron, fj, is the Fermi energy, defined 
as the average of first ionization and affinity energy. Param- 
eter values are U/t=8, V±/t=-6, J±/t=8, t±/t=l. The dot- 
ted line shows PES/IPES from the half-filled ground state 
(Stot~0), the full line from the half-filled 3 (ir, 7r)o state with 
S z =0 and S fot =l (cf. Fig. |) 
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FIG. 8. PES/IPES from the two-hole ground state 
°(0,0)_i (dotted line) and the half-filled 3 (7r,7r)o state with 
S z =S t ot= 1 (full line) of a 6- rung ladder (cf. Fig. |). The re- 
moved electron has .[-spin. The spectra for the half-filled state 
have been offset in y-direction so as to faciliate the compari- 
son. Parameter values are as in Fig. m. 



we have precisely the 4 bands expected from Fig. |5|. The 
physical interpretation is as follows: the initial 3 (tt, tt)q 
state carries a single triplet-like Boson in the "RVB- 
vacuum" . The "main bands" result from a creation of 
the photoholc in a singlet rung; then the final state, con- 
sisting of the initial triplet plus a propagating singly- 
occupied rung carrying the entire momentum transfer, 
belongs to the (2, 1) irreps. As the two excitations, triplet 
and hole, can scatter from each other, the main band be- 
comes broadened. But the photohole can also be created 
in the initial triplet, which results in a singly-occupied 
rung propagating in the RVB- vacuum, thus the new state 
belongs to the (1,0) irreps. Since in this second case 
the photohole has to absorb the momentum of the ini- 
tial triplet, (7r,7r), the sidebands' dispersion is shifted by 
(7r, 7r) with respect to the main bands. 

Proceeding to a 3 (ir, tt)q initial state with S z =l (in- 
stead of S* 2 =0) (Fig. |, full line), we get a PES/IPES 
spectrum similiar to the one for S z —0, but the sideband 
in PES has disappeared, whereas the one in IPES has 
gained some additional weight. The interpretation fol- 
lows from Fig. 0: As the initial state belongs to the (1,1) 
irreps, a PES transition into the (1,0) multiplet is now 
impossible. 

At this point, we remember that the SO(5) symmetry 
of our model implies that spin polarization and hole dop- 
ing are equivalent in that the empty rung is the "SO (5) 
partner" of the triplet rung. Consequently, Fig. |^ also 
shows the spectra originating from the two-hole ground 
state °(0, 0)_i (dotted line). Both states belong to the 
(1,1) irreps. Accordingly, their energies agree to com- 
puter accuracy. Since also the allowed final states of these 
two initial states belong to one and the same multiplet, 
namely (2,1), we expect their PES/IPES amplitudes to 
be directly related. A more detailed consideration ||] 
shows that their amplitudes are in fact identical - just as 
Fig. ^ demonstrates. The physical reason for this iden- 
tity is that the photohole cannot be created in the triplet 
rung, because the latter only contains two spin-up elec- 
trons. 

In summary our numerical "experiments" demon- 
strated, that SO (5) symmetry implies that the single par- 
ticle spectra in the doped ground state are identical to 
those of certain higher-spin states at half-filling. 



C. Recovering SO(5) symmetry at low energies 

The renormalization-group (RG) route has been 
proven to be particularly suited to describe the low- 
energy behavior of systems between one and two dimen- 
sions |34|-|36| . For this reason, we have carried out a RG 
study of a rather general ladder model of the Hubbard 
type The main result is that the effective Hamilto- 
nian of the system, i.e. the Hamiltonian which describes 
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excitations below a certain energy to*, becomes SO(5) in- 
variant for a wide range of bare models. These include 
interactions with quite arbitrary values of U, V, and V± 
(cf. Fig. ||), provided they are weak, as well as a moderate 
next-nearest-neighbor hopping t 2 . It is, thus, not neces- 
sary to introduce unphysical values for the parameters 
by hand in the model, since the exact SO (5) ladder dis- 
cussed in section [V.2 will be eventually recovered at low 
energies. This result is remarkable, since it shows that 
quite general Hubbard- like models, which are relevant for 
the description of the cuprate materials between ID and 
2D, although not explicitly SO (5) invariant on the micro- 
scopic level, display an exact SO (5) symmetry when ob- 
served on a macroscopic length scale oc Ef/uj* (in units 
of the lattice spacing) || (for the t 2 — case, see also 
Ref. ||). Physically, this means that the "standard de- 
viation" between the multiplet splitting (cf. Fig. ||) goes 
to zero in these models for low energies or large length 
scales. An important issue here is the introduction of a 
next-nearest-neighbor term t 2 , which is known to be es- 
sential in order to correctly describe AF correlations and 
Fermi-surface (FS) topology in cuprate materials. 

Specifically, we have considered two coupled chains 
with total low-energy Hamiltonian H = Hq + Hi, where 
Hq represents the non-interacting part and Hj the inter- 
action. This Hamiltonian describes interacting Fermions 
(expressed by creation and destruction operators Cko- and 
c k<r) c l° se to the FS. The FS consists of four points: two 
bands fcj^ = and fcj_ = tt, each one with two Fermi 
points corresponding to right- and left-moving Fermions. 
Since we are interested in low-energy properties very close 
to the FS, the Fermion dispersion can safely be taken as 
linear around the FS with (in general) band-dependent 
velocities vf- 

The idea of the RG is to divide the electronic excita- 
tions within the Brillouin zone into high-energy and low- 
energy modes, the latter being modes restricted within 
an energy lu <C Ep from the Fermi energy. One then 
eliminates the high-energy modes by integrating them 
out and constructs an effective Hamiltonian which is re- 
stricted to the low-energy excitations. The parameters 
of the new effective low-energy Hamiltonian thus depend 
on the energy cutoff u>. Also the total spectral weights at 
the FS (quasiparticle weights) Zk ± u> of the two bands are 
reduced due to the reduction of the cutoff and depend on 
ui. For this reason, in order to recover the canonical Fermi 
operators within the low-energy subspace, one should re- 
absorb Zk ±UJ into the definition of the Fermi operators 
and transform to c^ a = <J Zk_ L Jc\ s . a . This transforma- 
tion has to be done in order to preserve the sum rules 
for the total integrated spectral weight Z = 1 within the 
restricted subspace. The Ck CT now acquire the meaning 
of canonical operators with the correct anticommutation 
relations within the low-energy subspace. 

In practice, the integration of the high-energy modes is 



carried out by decreasing u> via infinitesimal steps, start- 
ing at to = Ep which corresponds to the bare (micro- 
scopic) Hamiltonian (see, e.g. p7j). Due to the restric- 
tion of the modes to a small "shell" around the FS, the 
scattering amplitudes in the interaction Hi can also be 
considered as dependent only on the Fermi momenta clos- 
est to where the corresponding processes take place as 
well as on the spin of the scattered particles (37|.|34|,|36| . 

The SO(5)invariant part Hj S °^ of Hi can conve- 
niently be written as a sum of products of SO (5) scalar 
contractions in the form ^>t vl/k, x where 
are SO(5) four-spinor Fermi operator (see Refs. [p8[|5|| 
for details). At w « Ep, the bare Hamiltonian is 
not SO (5) invariant in general and Hi consists of two 



parts Hi = H 



_ rr(SO(5)) 
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symmetry-breaking term. On the other hand, for t 2 = 0, 
the non-interacting part Hq is SO(5) symmetric, whereas 
a finite t 2 breaks PH (and thus SO (5)) symmetry also in 
Hq through a difference Avf between the Fermi veloci- 
ties of the two bands. 
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FIG. 9. RGnowoftheSO(5)-breakingtermsoftheHamil- 
tonian (in arbitrary units) ^j breakm s) (full), //( PH ~ brcakm s) 
(dashed), Avf (dotted) as functions of the RG flow 
r = — logui/ ' Ef = for U = 1., t = t± = 1, ti = 0.5 and half 
filling. The inset shows the t 2 = case where only //j broakm s) 
is nonvanishing. 

The striking result of our calculation is that, by inte- 
grating away the high-energy modes, from some energy 
scale u> = to* on, the SO(5)-breaking part of the Hamilto- 



nian ^jj brcakln s) eventually vanishes with respect to the 
SO(5)symmetric part. In the inset of Fig. 0, we show 
the RG flow for the t 2 = case (here the SO(5)-breaking 
term appears only in the interacting part of the Hamil- 
tonian). In the bulk of the same figure, the flow for the 
t 2 ^ case is displayed, for which there are three types 
of SO(5)-breaking terms, namely, Avf, a PH-breaking 
term _ffj PH blcakln s) anc [ a pi am SO(5)-breaking term 
^(breaking) ^ ^ three terms go to zero at to*. 

This result is of extreme importance, since it means 
that even though the system is not SO (5) invariant at 
the bare level, i. e. at ui = Ep, the SO (5) invariant 



11 



part of the Hamiltonian eventually dominates with re- 
spect to the symmetry-breaking part at energies smaller 
than u>*. We have checked that this occurs for very gen- 
eral values of the Hamiltonian, including on-site (U) and 
nearest- neighbor (V and V±) interactions, provided they 
are weak. 

An interesting result is that for the t% ^ case, the 
quasiparticle weights Zk ±UJ for the two bands renormal- 
ize differently. These quasiparticle weights have to be 
reabsorbed into the definition of the canonical Fermi op- 
erators through the transformation to the Cko- variables, 
as explained above. This has the important consequence 
that the new low-energy SO(5) invariant Hamiltonian 
jj(SO(5)) j g invariant under a renormalized SO(5) 
symmetry in terms of new fr operators (cf. (0) ) , whereby 
the Fermi operators Cko- are replaced with Ck<r- This ex- 
tended concept of SO (5) symmetry makes it possible for 
this generalized symmetry to occur in a larger and more 
generic class of physical systems than the ordinary SO (5). 
Moreover, a SO(5) theory unifying antiferromagnetism 
and superconductivity in terms of this generalized repre- 
sentation can admit possible asymmetries between the 
antiferromagnetic and superconducting phase, like for 
example the difference in T c or in the order parameter 

In conclusion, the renormalization-group study shows 
that the effective low-energy Hamiltonian of a quite 
generic Hubbard-like ladder with weak interaction is 
SO (5) symmetric. This holds true also with the inclu- 
sion of a next-nearest-neighbor hopping t%, provided it 
is written in terms of the appropriate canonical Fermi 
operators c for the low-energy subspace. 
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